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Abstract: This paper presents an analytical study for the magnetohydrodynamic (MHD) flow of a 
generalized Burgers 'fluid in an annular pipe. Closed from solutions for velocity is obtained by using finite 
Hankel transform and discrete Laplace transform of the sequential fractional derivatives. Finally, the 
figures are plotted to show the effects of different parameters on the velocity profile. 
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I. Introduction 

Fluid non-Newtonian does not show a linear relationship between stress and strain rate and received a 
lot of attention due to the increasing applications of industrial and technological field, such as polymer 
solutions, paints, blood and heavy oils. No model alone that can describe the behavior of all fluids non- 
Newtonian because of complex behavior. 

Thus, it was suggested many of the foundational equations of non-Newtonian fluid models. The 
Burgers' fluid is one of them which cannot be described a typical relation between shear stress and the rate of 
strain, for this reason many models of constitutive equations have been proposed for these fluids[7,8,ll]. Many 
applications of this type of fluid can be found in [1,2,15,18]. Many of the developments in the theory of 
viscoelastic flows have been mainly restricted to the formulations of the basics equation and constitutive models 
[12,16], and many applications of fractional calculus can be found in turbulence and fluid dynamics, stochastic 
dynamical system and nonlinear control theory [3,14,19]. 

Recently, Tong and Liu [6] studied the unsteady rotating flows of a non- Newtonian fluid in an annular 
pipe with Oldroyd- B fluid model. Tong ... etc [5] discussed the flow of Oldroyd- B fluid with fractional 
derivative in an annular pipe. Hyder ... etc [17] discussed the flow of a viscoelastic fluid with fractional 
Burgers' model in an annular pipe. Tong ... etc [4] discussed the flow of generalized Burgers' fluid in an 
annular pipe. Khan [13] investigated the (MHD) flow of generalized Oldroyd- B fluid in a circular pipe. Later 
on [10] investigated the slip effects on MHD flow of a generalized Oldroyd- B fluid with fractional derivative. 

In this paper, our aim is to steady the effects of MHD on the unsteady flow of a viscoelastic fluid with 
fractional generalized Burgers' fluid model in an annular pipe. The exact solution for velocity distribution is 
established by using the finite Hankel transform and discrete Laplace transform of the sequential fractional 
derivatives. 

II. Governing Equations 

The constitutive equations for an incompressible fractional Burger's fluid given by 

t = pi + s , (i + rfif + rb 2 t a >s = + ^5f >a (i) 

where T denoted the cauchy stress, -pi is the indeterminate spherical stress, Sis the extra stress tensor, 

A = L + L is the first Rivlin- Ericksen tensor with the velocity gradient where L = grad V , ju is the dynamic 
viscosity of the fluid, A l and X 3 (< X x ) are the relaxation and retardation times, respectively, X 2 is the new 
material parameter of Burger's fluid, a and ft the fractional calculus parameters such that 0 < a < j3 < 1 and 

Df the upper convected fractional derivative define by 

D?S = D?S + (V.V)S - LS - SL T (2) 

DfA = DfA + (V.V)A LA AL 1 (3) 

in which D? and Df are the fractional differentiation operators of order a and ft based on the Riemann- 
Liouville definition, defined as 
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D p [f (01 = ^ } dr ,0 < p < 1 (4) 

r(i-p) a j (t-ry 

here T(.) denotes the Gamma function and 

D 2p S = D p (D p S) (5) 

The model reduced to the generalized Oldroyd- B model when A 2 = 0 and if, in addition to that, 

a = J3 = 1 the ordinary Oldroyd- B model will be obtained. 

We consider the MHD flow of an incompressible generalized Burger's fluid due to an infinite accelerating 
plate. For unidirectional flow, we assume that the velocity field and shear stress of the form 

\ = w(r,t)e z ,S = S(r,0 (6) 

where e z is the unit vector along z - direction .Substituting equation (6) into (1) and taking account of the initial 
condition 

S(r,0) = 0 (7) 
we obtain 

+ +r 2 D 2 t a )S vz = //(1 + vlfDf )3 r *KM) (8) 

(1 + W +^D t 2a )S zz -2S a (% +^Df )S r w(r,0 = -2^(d r u<r,0) 2 (9) 
where =S 0Z =S r0 =S 00 = 0. Furthermore, it assumes that the conducting fluid is permeated by an imposed 
magnetic field H = [0, H 0 , 0] which acts in the positive 6 - direction. In the low- magnetic Reynolds number 



approximation, the magnetic body force is represented as ctHqW , where <j is the electrical conductivity of the 

fluid. Then in the presence of a pressure gradient in the z- direction, the equation of motion yields the following 
scalar equation: 

dw dp 1 d tt2 

at oz r or 

where p is the constant density of the fluid. Eliminating between Eqs. (8) and (10), we obtain the following 
fractional differential equation 

(1 + ylfD? + X^l\ a ) — = - - (1 + ^ + %p\ a ) ^ + 

a p 2 dz (11) 

rsZ 1 Pi 

+ v(l + ^Df X-^y + - —)w - M(l + 4 a D? + %D* a )w 

dr r dr 

crH 2 

where v = — is the kinematic viscosity and M = — is the magnetic dimensionless number. 

P P 

III. Plane Poiseuille Flow 

Consider that the flow problem of an incompressible generalized Burgers' fluid is initially at rest 

between two infinitely long coaxial cylinders of radii R 0 and ^ ( > R 0 ). At time t = 0 + the fluid is generated 

due to a constant pressure gradient that acts on the liquid in the z- direction . Referring to Eq. (11), the 
corresponding fractional partial differential equation that describe such flow takes the form 



+ +r 2 D 2 t a ) — = -A 

dt 



+ 



/ , -a .-2a \ 

i+r — — +r 2 — 

v r(i-«) r(i-2a) y 

rs2 IP) 

+ v(l + AfDf)(-^- + -— )w-M(l + AfDf +r 2 D 2a )w 

dr r dr 



(12) 



where A = — — is the constant pressure gradient 
p dx 



The associated initial and boundary conditions are as follows 

^w(r,0) = ^ 
dt dt 

w(R 0 , t) = w(R { ,0 = 0 ,J>0 (14) 



w(r,0) = — w(r,0) = — w(r,0) = 0 , R 0 < r < R l (13) 



To obtain the exact analytical solution of the above problem (12)- (14), we first apply Laplace 
transform principle [9] with respect to t, we get 
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S (i+A^ a +r 2 s 2a )w = -— (\+r lS a +r 2 s 2a )+ 



d 2 | 1 S„ 
dr 2 r dr 



+v(i+a^x-_-t+-— )w-M(i+r lS a +r 2 s 2a )w 



(15) 



w(r,0) = 0 

vv(7? 0 , 5-) = vv(7?, , s) = 0 , f > 0 

where w(r, s) is the image function of w(r, t) and s is a transform parameter. 
We use the finite Hankel transform [9], defined as follows 



(16) 



w H = j r w B 0 (rk t )dr 



I • • • 



(17) 



*0 

and its inverse is 

2 oo 7_2 — 



_ 71 
W = 



Z 



(18) 



2 / 0 2 (*o*/Wo(*i*/) 
where k { are the positive roots of equation B 0 (R 1 k i ) = 0 and B 0 (rk t ) = J 0 (rk i ) Y Q (R 0 k i ) — Y 0 (rk t ) J 0 (R 0 k i ) 
where J 0 (.) and Y 0 (.) are the Bessel functions of the first and second kinds of order zero. 
Now applying finite Hankel transform to Eqs. (15)-(16) with respect to r, we get 

A (1 + A^s a +A a 2 s 2a ) 



w H = 



s (s + M)(1 + a^ +r 2 s 2a ) + v 1 k 2 {l + ^s p ) 



(19) 



Now, writing Eq. (19) in series form as 



oo 



a,b,c,d ,n>0 



W 



H 



= - A(l + 5 a + A a 2 s 2a (- 1) * ^k\M c+d (v k?) n (A") 



-k-l+d 



k=0 



a+b+c+d+n=k 



s 



n 



(20) 



f ..t_2 ,,\ k+l 



a\b\c\d\n\ 



a+l V k i M 

S + — - + — 



V 



J 



na 

where S = k — 1 + 2a(b + c)-{c + d + n)- aid + n) + fin-an . And its discrete inverse Laplace transform [20] 
will take the form 



X 



a,b,c,d,n>0 



Mc+d / . 7 2\n / na \ —k—l+d / o# \£>+c / n0\n 
(yki ) (\) (/I,) (4) . (a+ i) t+ ( a+ i-5)-i 



w a =-Aj;(-i)' x 



' a+\,a+\-d 



V 



"^ (l *' 2+M)f 



a+l 



5 or 



a\b\c\d\n\ 

f 



J 



a+\,a+\-8 



V 



a+l 



2a a+l,a+l—S 



J 



V 



"^ (l *' 2+M)f 



a+l 



J 



(21) 



oo 



where E^(z) = ^- 



(j + m)\z\ 



j\T(aj + am + P) 
following property of inverse Laplace transform is used [20] 



is the generalized Mittag- Leffler function [9] and to obtain Eq. (21), the 



L < 



mis 



(s A Tc) 



y = t Am+M ~ 1 E™ M (±ct A ) 



f 



i A 



Re(s) > 



V 



(22) 



J 



Finally, the inverse finite Hankel transform gives the analytic solution of velocity distribution 



w(r, t) = - 



A* 2 ^ k? B 0 (rkMR&) 



2 j-fJli^-JliR^) 



00 



I*" 1 )* I 



a,b,c,d ,n>0 H/fC+d , i 2\n / na\-k-l+d / na\b+c/ n/3\n 

M {vk t ) {\) (/{^ (aJ) Aa+1)kHa+1 -S)-l 



k=0 



( 



'a+l,a+l-S 



V 



a+l 



J 



+ ^E k 

^a a+l,a+\-8 



a+b+c+d+n=k 
( 



a\b\c\d\n\ 



f 



V 



a+l 



+ ^E k 



f 



J 



t 



2a a+l,a+l-S 



V 



a+l 



J 



(23) 
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3.1 The limiting cases 

1- Making the limit of Eq.(23) when a ^ 0 and M — > 0 (c=d=0) , we can get similar solution velocity 
distribution for unsteady flows of a viscelastic fluid with the fractional Burgers' model, as obtained in Ref[17]. 
Thus the velocity field reduces to 



w(r t)= A * 2 f k > BMiVlm 

2 £/ 0 2 (/^Wo(W 



oo 



J a+l,a+l-S 



V 



1 



a 



J 



+ ^-E k 

a+\,a+\-S 



E(-i) 1 I 

k=0 a+b+n=k 
f 



a,b,n>0 / ]A\ n ( ; a / Q « \b-n / ; /? \n 

(V/c. ) ) (A 2 ) _{A 2 J_ t(a+l)kHa+l _ S) _ l 



1 



V 



alblnl 



1 a 



J 



t 



2a a+l,a+\-S 



( 



V 



1 



J 



(24) 



where 5 = k-l + 2ab-n-2an + /3n. 



2- Making the limit of Eq.(23) when a ^ 0 , >l 2 — » 0 (b=0) and M — » 0 (c=d=0) , we can get the velocity 
distribution for a generalized Oldroyd- B fluid. Thus the velocity field reduces to 



w(r, t) = - 



Att 2 ^ kf B^rk^JliR^) 



2 / 0 W,)-WA) 



a ,n>o (vk 2 Y ( 2 a \~ k ~ l ( A^Y 



£=0 



V 



2 



+ 



a+n=k 
f 



a 



' a+l,a+l-S 



J 



V 



2 \ + a+l 



J 



(25) 



where S = k-l-n + /3n . 



IV. Numerical Results And Discussion: 

In this work, we have discussed the MHD flow of generalized Burger's fluid in an annular pipe. The 
exact solution for the velocity field U is obtained by using the discrete Laplace and finite Hankel transforms. 
Moreover, some figures are plotted to show the behavior of various parameters involved in the expressions of 
velocity U . 

A comparison between the magnetic parameter effect M (Panel a) and no magnetic parameter effect 
(M=0) (Panel b) is also made graphically in Figs 1-5. 

Figs. 1 and 2 provide the graphically illustrations for the effects of the non- integer fractional parameter 
a and (3 on the velocity fields. The velocity is increasing with the increased the a and /? for both cases 

(M = 0&M^0). 

Fig. 3 provides the graphical illustration for the effect of relaxation parameter \ on the velocity fields. 
The velocity is decreased with the increase of \ for both cases (M = 0& M^0). 

Figs. 4 and 5 are prepared to show the effect of the material parameter A 2 and the retardation parameter 
/^on the velocity field. The velocity is increasing with the increase of A 2 and X$ for the both cases 
(M = 0&M^0). 

Fig. 6 is established to show the behavior of the magnetic parameter M for small as well as for long 
time. It is observed in (Panel a) that for short time t = 0. 1 the increase in magnetic field M will decrease the 
velocity profile, while quite the opposite effect is observed for long time t = 0.5 in (Panel b) i.e., the increase in 
magnetic field M will increase the velocity profile. 

Comparison shows that the velocity profile with magnetic field effect is larger when compared with the 
velocity profile without magnetic field effect. The effect is explained on the long time. 
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Fig. 1. The velocity for different value of a when keeping other parameters fixed a) M = 3 b) M = 0 




Fig. 2. The velocity for different value of p when keeping other parameters fixed a) M = 3 b) M = 0 




1.2 1.4 1.6 1.8 2.0 1.2 1 4 1.6 1.8 2.0 



Fig. 3. The velocity for different value of \ when keeping other parameters fixed a) M = 3 b) M = 0 
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Fig. 4. The velocity for different value of X 2 when keeping other parameters fixed a) M = 3 b) M = 0 




Fig. 5. The velocity for different value of A 3 when keeping other parameters fixed a) M = 3 b) M = 0 




Fig. 6. The velocity for different value of M when keeping other parameters fixed a) t = 0.1 b) t = 0.5 
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